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Summary. We review the basic ideas lying at the foundation of the recently 
developed theory of twisted symmetries of differential equations, and some of 
its developments. 

Introduction 

The study of nonlinear differential equations was the main motivation to Sophus 
Lie when he created what is nowadays known as the theory of Lie groups and 
Lie algebras, which is by now recognized as one of the most effective tools for 
the study of differential equations - both in geometrical sense and for what 
concerns the search for explicit solutions [25, 37, 55, 59]. 

The attention of Lie was mainly focused on what we now call Lie-point 
symmetries, but very soon generalizations of these were investigated, starting 
with contact symmetries. The effectiveness of symmetry methods for differen- 
tial equations - greatly increased in recent years by the possibility to perform 
the often involved computations required by these via symbolic manipulation 
languages, i.e. by computer - led various authors to consider generalizations in 
several directions [15, 16, 23]; all of these however - at the exception of one, i.e. 
the subject of this paper - are based on the same scheme: one considers a vector 
field X acting on the full phase manifold M - the manifold of independent and 
dependent variables 1 - and then takes into account its action on the jet man- 
ifold J k M of suitable order (i.e. with k the order of the differential equation 
one wishes to study), which is obtained by the standard prolongation procedure. 
That is, once we know how X acts on the independent variables x % and on the 
dependent ones u a , we also know how it acts on the derivatives (of any order) 



1 In most symmetry applications, this is just a linear space M = B X U with B and U 
the spaces of independent and dependent variables respectively. However equations can be 
defined on a manifold, in which case B is not a linear space, or dependent variables take 
value in a nonlinear manifold U, so that M = B X U is also not such; or the equation can be 
complemented with side conditions, in which case M is in general a bundle over B. 
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of the u's with respect to the x's. In computational terms, this is readily de- 
scribed by the usual prolongation formula [25, 37, 55, 59]; in geometrical terms, 
this makes use of the contact structure which naturally equips the jet manifold 
J k M [8, 56]. Thus, the generalizations lie in the class of admitted vector fields 
X on M, i.e. on what qualifies (if satisfying certain conditions specific to the 
equation under study) as a "symmetry" ; while once this is given, the action on 
higher order jet bundles J k M is just the natural prolongation of the one in M; 
in other words, we have the vector field X^ - the natural prolongation of X - 
acting in J k M. 

Here we are interested in the exception mentioned above, i.e. in what we 
would like to call twisted symmetries (the reason for this name will be clear 
in the following). These were first proposed by M.C. Muriel and J.L. Romero 
under the name of "C^M^ )-symmetries" or "C°°-symmetries" for short, or 
finally of "A-symmetries" , in the framework of ODEs [41, 45, 46, 47, 57], and 
then generalized to PDEs under the name of "^-symmetries" [22, 26, 29]. 

It should be stressed that in this case the generalization with respect to 
standard Lie symmetries lies not in the definition of symmetry 2 , but in the 
way vector fields are prolonged from the phase bundle M to jet bundles J k M. 
In this, they arc deeply different from other proposed generalizations of Lie 
symmetries. 

Definition and properties of twisted symmetries will be discussed below, 
but we like to recall immediately that these allow to integrate - by symmetry 
methods - equations which are known not to admit any standard Lie symmetry 
[41, 43, 44]; needless to say, this is a major reason (but not the only one) for the 
interest they raised. It also turned out that twisted symmetries are interesting 
from the point of view of the geometry of differential equations. Here we will try 
to keep some balance between these two points of view, the (mostly analytic) 
applicative and the (mostly geometric) theoretical one. 
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Cicogna and P. Morando, who shared my investigation on twisted symmetries; 
and D. Catalano, M.E. Fels, F. Magri, G. Marmo, P.J. Olvcr, G. Saccomandi 
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is a pleasure to thank the Editors for their invitation to contribute a paper to 
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2 I.e. in the class of acceptable vector fields or in the equation-related specific conditions 
to be satisfied by the prolonged vector field. 
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1 Standard prolongations and symmetries 



In this section we will very briefly recall the standard (and well known) notions 
of prolongation and symmetry, to the aim of fixing notation. 

We suppose the reader has some familiarity with symmetry of differential 
equations and hence do not enter in details; these can be obtained e.g. from 
[25, 37, 55, 59]. Summation over repeated indices will always be assumed unless 
otherwise explicitly stated. 

1.1 Prolongations 

As already stated, we will denote the phase bundle as M = B x U, where B 
is the manifold to which the independent variables belong, equipped with local 
coordinates x l (i — l,...,q); and U is the manifold in which the dependent 
variables (or fields) take values, which is equipped with local coordinates u a 
(a= l,...p). 

We consider a (Lie-point) vector field X in M; this will be written in coor- 
dinates as 

x = ft'^w + ^ u) i- (1) 

As well known, given a section oo = (x,u = fo(x)) of the bundle M — B xU, 
this is transformed under the infinitesimal action e eX of X into the section 
o e = {x,u = f e {x)) with 

f e (x) = f (x) + S (p-r<)0 , 

where the subscript "0" reminds that the functions ip and £ - as well as the 
derivatives uf - should be computed on cfq, i.e. 

p°L = v a (x,f (x)) , Clo = C(x,fo(x)) , uXo = (dfo(x)/dx) . 

This is also said by considering the evolutionary representative (or vertical rep- 
resentative) Xy Of X, 

x v = Q a ^ , withg Q - <p a (x,u) - KC(x,u) . 

The jet bundles J k M associated to M can be seen also as bundles over B, 
J k M = B x (here is the manifold of fields and their derivatives of 
order up to k; we will denote by the manifold of field derivatives of order 
exactly k). They can also be seen as bundles over the jet bundles of lower order, 
J k M = J fe_1 M x U^ k \ and it is this structure which comes into play when 
considering the prolongation operation in recursive terms. 

The first jet bundle JM is naturally equipped with local coordinates (x l , u a , uf). 
If X is given in coordinates by (1), its first prolongation is 
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where the coefficients ipf are given by (see below for the derivation of this 
formula) 

C = D^ a - u a 3 D£ . (3) 

Similarly, the jet bundles J k M are equipped with local coordinates (x l , u a , Uj) 
with J multi- indices of order up to k. If X is given in coordinates by (1), its 
fc-th prolongation is 

Xik) = 5? + ^ U ^a + E rjM^, (4) 

IJI=1 



where the sum is over all multi-indices J of order 1 < | J < k, and the coefficients 
ipj are given by the recursive prolongation formula (here one should understand 

ij>% = D^j - u a Jtk Die . (5) 

Note that, as well known, for vector fields in evolutionary (vertical) form, 
the prolongation formula reduces to 

r Jti = ii'-:, . 

In the case of ODEs, when we have only one independent variable x, the 
above notation is a bit too heavy, so it is convenient to also have a simpler 
notation for this case. We will denote the derivative of order j as 

u « := (d 3 u a /dx 3 ) , 

and the fc-th prolongation of a vector field X = £(x, u)(d/dx) + ip a (x, u)(d/du a ) 
will be written as 

with the V'j satisfying 

V>? +1 = £> X V? - u a J+1 D x i . (7) 

Last but not least, we recall that given any two vector fields in M, the 
commutator of their prolongations is the prolongation of their commutator, i.e. 

XW,YW] = ([X,Y]) (k) . (8) 



The prolongation formula analytic derivation 

We would like to recall how the prolongation formula is obtained in analyt- 
ical terms; in a later subsection we will discuss its geometrical meaning and 
correspondingly a geometrical derivation. 
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Let us consider a section a = (x,u), corresponding locally to a function 
u = f(x); let x\ and x 2 = x\ + 8x % be two nearby points in B, with Sx l an 
infinitesimal displacement in the direction x % . The corresponding points on a 
will be pi = (x\,ui — ,f{x)) and p 2 — {x 2 ,u 2 = /(a^)); the partial derivatives 
of u on a are given by 



ut{x) = lim := lim R. 

Sx^O X 2 - X\ 

Needless to say, u 2 = + u^8x l + o{8x l ). 

If now we act in M by X, i.e. by its infinitesimal action e eX , the points pi 
and P2 are mapped to new points p\ = (xi,ui) and p 2 = (x 2 ,u 2 ) with 



x j ~ x j + e £ ( x j> u i) 



Uj + Eip(Xj,Uj) . 



We have now to compute (the limit of) 

- _ u 2 - ui [u 2 + eifi(x2,u 2 )] - [ui + e{x 1 ,ui)] 

X2 - xi [x 2 + e£(xi,ui)] - [xi + £(xi,Ui)} 

(u 2 - ui) + e[ip{x 2 ,u 2 ) - (p(xi,ui)] 
(x 2 - x\) + e[£_(x 2 ,u 2 ) - £Oi,«i)] 

It results immediately, with all derivatives evaluated in (x\,y\), 



C(x 2 ,u 2 ) = r(zi,ui) + ^fe l + ^<&r' «i) + [Di(Z)]6x l ; 



ip a (x 2 ,u 2 ) 



dx 

dtp a ■ 
tp a (x uUl ) + -ttjSx 1 
ox 1 



du b 

+ -^u°5x l = ip(x!, Ul ) + [Di(ip)} 8x l 



Here and below, Di denotes the total derivative with respect to x 

_d_ 

dx 1 



D, 



We have therefore, recalling (x 2 — x\) = 8x l , writing (u 2 
short, and omitting terms of higher order in e, 



(9) 

Mi) = Su for 



R 



8u + e[(D t Lp)8x l 



Su + e[(Di(p)5a 



{x 2 - xt) + e[(AO^*] 5x% I 1 + e(A01 
\8u + e\(Di<p)8x*] [l-g(AQ] ^ g 

<5x* 



i? + £ 



(Av) - ^ (AO 



(5u 



r(AO 



This yields, upon taking the limit <5x l — >• 0, precisely the prolongation formula: 
if X is given in coordinates by (1), its first prolongation is 



du° 



(10) 
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where the coefficient ipf is given by 

C = D^ a - u a 3 D^ . (11) 

The computation for higher order prolongation is exactly the same (with the 
role of u a taken by u"). 

1.2 Symmetries 

As well known, a differential equation E of order k defines a submanifold Se 
(of codimension one) in J k M, also called the solution manifold for E. The 
same applies for a system A = {E 1 , ■■■ } E P } of equations; in this case Sa is of 
codimension p (if the system is non degenerate). 

If the vector field X on M is such that its fc-th prolongation X^ leaves 5a 
invariant, i.e. 

X™ : S A -> TS A , (12) 

we say that X is a (Lie-point) symmetry of A. The set of symmetries of A will 
be denoted as Xa- 3 

It follows immediately from this definition - and from (8) - that the sym- 
metries of a given A are a Lie algebra (under the commutator): 

X G Xa , Y G Xa => [X, Y] G X A . (13) 

In order to check if a given vector field X is a symmetry of A, one just has 
to check that 



X {k \E a ) 



= . 

Sa 



1.3 The geometric meaning of prolongation 

The prolongation operation was defined before as the lifting of the action of a 
vector field on independent variables and fields, to the derivatives of the latter 
with respect to the former. In this way, it was defined in an analytic way. 
But, prolongation has a deep and intrinsic geometric meaning as well, and it is 
convenient to focus on this in order to consider, later on, twisted prolongations. 

The first jet bundle JM is naturally equipped with local coordinates (x z , u a , uf); 
in terms of these, the contact structure 9 1 is described by the contact forms 

6 a = du a - uf dx* . (14) 

The prolongation X^ of X is the unique vector field in J 1 M which coincides 
with X on M and preserves the contact structure 1 . 

Similarly, the jet bundles J k M are equipped with a contact structure 6 fc ; in 
terms of the local coordinates (x l ,u a , Uj), the contact structure <d k is described 
by the contact forms 

6 a j = du a j - u'jjdx* , (15) 

3 More precisely, X is a symmetry generator, and e aX is a one-parameter symmetry group 
for A. This slight abuse of notation is commonplace in the literature, and we will keep to it. 
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where < \ J\ < k — 1. The prolongation X^ of X is the unique vector field in 
J k M which coincides with X on M and preserves the contact structure fe . 

It is maybe worth defining more precisely, for the benefit of the reader less 
used to this geometric language, what is meant by "preservation of the contact 
structure" . 

Consider a set $ r } of generators for fc , and denote by £ the C°°(J k M) 

module generated by these same generators 4 . We say that a vector field Y on 
J k M preserves 8 fe if and only if 

C Y (ti)££ Wg£; (16) 

here and in the following C Y is the Lie derivative under Y. 

Note, for later reference, that the condition (16) can be expressed equiv- 
alently in terms of conditions involving the commutator of Y with the total 
derivative operators in particular, it is equivalent to either one of 

[Di,Y]_\0 = W G £ , (17) 
[Di,Y] = h?D m +V , (18) 

with h™ G C°°{J k M) and V a vertical vector field in J k M (seen as a bundle 
over J k ~ 1 M, see above). 



The prolongation formula geometric derivation 

We have recalled above how the prolongation formula is obtained in analytic 
terms; here we want to show for comparison how it is derived in geometric 
terms, i.e. from (16). We will again limit to consider the first prolongation, 
computations being the same for higher ones. 

First of all we note that, from the general properties of the Lie derivative, we 
have just to show that C Y {6) G £ for all 9 in the generating set of the contact 
structure. In facts, for d = /3 a 6 a with j3 a G C°°(JM), we have 

c Y (M a ) - [c Y {PaW + p a [c Y {e a )] ■ 

the first term in the right hand side is by definition in £, and the second is 
always in £ if and only if £y(# a ) G £ for all a = 1, ...,p. 
In our case, d6 a = dw? A dx l . We will write 



so that Y_l9 a = ip a - ufC and d(r_l6» a ) = dip a - <df - f duf; moreover 
Y_ld0 = Cduf-tpfdx\ Therefore, using the Cartan formula C Y (6) = r_ld6»+ 
d(Y~_l#), we get immediately 

£ Y (8 a ) = dip a - u^dC 3 - ip?dx l . 



4 I.e. the set of one-forms obtained as sum of the i?j with coefficients which are C°° functions 
on J k M, 8 = YZ=l ^(x,^^)^. 
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Using du a = 9 a + u*dx\ we have dip a = (d(p a /du b )0 b + (D^ a )dx 4 and d£» = 
(d£ j /du b )6 b + (D i ^)dx i ; thus we can rewrite 

£ Y (6 a ) = [(d<p a /du b ) + (de/du b )]6 b + [(D^ a - v%Di?) - V?] dx l . 

The first term on the right hand side is surely in £, by definition, while the 
second is either zero or surely not in £. We conclude that £y(6 |a ) E £ if and 
only if 

C - (D^) - u« (A^) • 
We have thus obtained again the prolongation formula (2). 

2 Lambda-prolongations and symmetries 

The generalization of classical Lie-point symmetries we wish to consider, i.e. 
twisted symmetries, modifies the standard prolongation operation. That is, one 
considers vector fields Y in J n M which are not the prolongation of some vector 
field X in M, i.e. Y ^ X^> for any X , yet are related to such a vector field 
in a precise manner to be discussed below. Quite surprisingly, when these are 
symmetries of a differential equation A - i.e. when Y : Sa — > TSa - they are 
still effective in obtaining symmetry reductions (for ODEs) or invariant solutions 
(for PDEs) of differential equations. 

We will start by discussing the first class of twisted symmetries to be dis- 
covered, i.e. A-symmetries. 

2.1 The work of Muriel and Romero 

In 2001, C. Muriel and J.L. Romero [41], analyzing the case where A is a scalar 
ODE, noticed a rather puzzling fact. 

They substitute the standard prolongation formula (7) with a "lambda- 
prolongation" formula 

= (D x + A)* fc - u k +i (D x + A) £ ; (20) 

here A is a real C°° function defined on J X M (or on J k M if one is ready to deal 
with generalized vector fields). For A = one recovers standard prolongations. 

We say that X is a "lambda-symmetry" of A if its "lambda-prolongation" 
Y is tangent to the solution manifold Sa C J n M. 

Then, as mentioned above, it turns out that "lambda-symmetries" are as 
good as standard symmetries for what concerns symmetry reduction of the 
differential equation A and hence determination of its explicit solutions. As 
pointed out by Muriel and Romero, it is quite possible to have equations which 
have no standard symmetries, but possess lambda-symmetries and can therefore 
be integrated by means of their approach; see their works [41, 43] for examples. 

It is quite remarkable that the A-symmetries approach is able to explain 
the - rather puzzling - fact that there were equations explicitly integrable by 
quadratures and not possessing any (standard) symmetry [4, 31, 32, 55]. 
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2.2 The invariants-by-differentiation property 

The possibility of using A-symmetries as effectively as standard ones for order 
reduction of differential equations has its origin in the fact the recursion formula 
allowing to build higher order differential invariants from lower order ones [55, 
56] , also called "invariants-by-differcntiation" , applies for A-symmetries as well. 
This point of view was stressed by Muriel and Romero in another paper [47], 
where other examples are also provided. 

If rj and £ are differential invariants for the vector field Y, the invariants-by- 
differentiation property states that p := (D x ()/(D x r]) (with D x r] ^ 0) will also 
be a differential invariant. The proof of this property goes as follows: acting 
with Y on p we have 

Y (2A\ _ [Y(D x Q}(D x y) - (D x Q[Y(D x y)} 
\D xV ) (D xV y 

Thus p is a differential invariant if and only if 

[Y{D x Q]{D x r,) = (D x Q[Y(D x r])] ; 



we can rewrite this equation as 

(D x r,)(\Y,D x ](Q)-(D x (YQ)(D x ri) 

We assumed r] and ( are differential 
and the above equation reduces to 

(D xV ) ([Y,D x }(0) 



= (D x C)([Y,D x }( V ))-(D x 0(D x (Y V )) . 
invariants for Y, hence Y(i]) = Y(() = 0, 

- (d x () ([r, £>*](»?)) • (2i) 



If Y is a A-prolongation, it satisfies (up to a term which vanishes on the contact 
distribution) 

[D X ,Y] = XY + hD x (22) 

with h a smooth function; with this, and recalling again Y(rf) — Y(Q — 0, the 
above equation reduces to the trivial identity 

h(D xV )(D x Q = (D x £)h(D x ri) . 
2.3 Systems of ODE 

Muriel and Romero also considered how the concept of A-symmetry would ex- 
tend to systems of ODEs [45] . 

They considered a A-prolongation given by 

r k+ i - (D x + XW k - 4+i (Ac + A)£ , (23) 

i.e. the immediate generalization of (20), and studied a system A of first order 
ODEs (which we will also call a dynamical system), 

u% = F a (x,u\...,u r ) , (a = l,..,r) . 
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In this case, if X is a A-symmetry of the system A, there exists a change of 
variables (y, w) = (3(x, u) under which the system reduces to a system of (r — 1) 
first order ODEs, 

wl = G?{x,w s ;-,w r - x ) , (a = l,...,r-l) 

and the variable w r must satisfy an auxiliary ODE 

H(y,w r ,w r y ) = 0. 

2.4 Lambda symmetries and nonlocal standard symme- 
tries 

In their seminal work, Muriel and Romero also remarked (see sect. 5 of [41]) that 
there is an intriguing relation between A-symmetries and exponential symmetries 
of differential equations. 

Exponential symmetries [55] represent a specific form of nonlocal vector 
fields, 

x = Jp(x,u)dx (^( x ,u)-^ + ^ := e f P{x ' u)dx X . (24) 

Needless to say, their prolongation X* must satisfy the usual symmetry condi- 
tion X* : S A TS S . 

Muriel and Romero proved that if X given by (24) is an exponential sym- 
metry for A, then X defined in (24) is a A-symmetry for 8, with 

A = P(x,u) . 

As nonlocal symmetries can be used, pretty much as standard symmetries, 
to reduce and integrate differential equations [7], this shows that all the reduc- 
tion/integration methods based on nonlocal symmetries (of exponential type) 
[5, 6, 7, 55] can be automatically formulated in terms of A-symmetries. 

In a later work [48], Muriel and Romero went further in their study; in 
particular they studied the connection between A-symmetries and the so called 
"type I hidden symmetries" (see below). The key property of use here, shown in 
their work, is that if an equation d™w/dx™ = f(x, u^ n ~^) admits a A-symmetry 
X, then when the equation is written in new variables (y, v) as 

d%/dy" = g{y,v {n ' 1] ) , 
the vector field X is a A-symmetry for the equation, with 

a * 

D x y 

Muriel and Romero focused in particular on first order equations; in this case 
a theorem by Adam and Mahomed [6] gives a criterion which must be satisfied by 
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a nonlocal symmetry of a first-order equation in order that a symmetry-related 
(computable) transformation maps the latter into an integrable equation; this 
theorem can be reformulated (in a simpler way) in the language of A-symmetries 
[48] . This result was then applied to relevant equations such as Riccati equations 
and Abel equations of the second kind [48]. 

The relation between A-symmetries and nonlocal standard symmetries was 
also investigated by Catalano-Fcrraioli [12] by means of the theory of coverings 
[37, 36, 61]. In this, one embeds J k M into a space of higher dimension by adding 
one or more extra variables w a with their (first) derivatives, and augments the 
equation under study with extra equations 

dw a /dx = h a {x,u {n \w) ; 

nonlocal symmetries can then be expressed as local symmetries of the augmented 
system. The connection with A-symmetries arises (with a single extra variable 
w) when one chooses h = A. We refer to the original paper [12] for details. 

2.5 Lambda symmetries and other types of symmetries 

Lambda symmetries turned out to be relevant, and providing a way to soundly 
understand, also other kinds of symmetries. Here we will just briefly mention 
these interrelations, referring the reader to the original papers. 

Gandarias, Medina and Muriel [30] discussed the relation between A-symmetries 
and the so called potential symmetries of differential equations [11]. They also 
discussed how this relation can be of help in the integration of differential equa- 
tions not possessing Lie-point symmetries. The same problem was also tackled 
in a paper by Muriel and Romero [46] , where they also discuss the relation with 
potential symmetries. 

Similarly, A-symmetries are instrumental in integrating equations with non- 
solvable symmetry algebras; this point was discussed by Muriel and Romero 
both in general terms [43] and for specific algebras relevant in Physics (and in 
other applications as well) [42, 44]. 

Lambda symmetries also bear some interesting relation to so called hidden 
symmetries (such as symmetries which got lost in the reduction process) [1]. 
See [2, 3] for a discussion. 

Hidden symmetries also have relations [33] with so called solvable struc- 
tures [9, 10]; these in turn have been very recently studied in connection to 
A-symmetries [13]. We refer the reader to the original paper [13] for this mat- 
ter. 

2.6 Lambda symmetries and integrating factors 

It turns out that A-symmetries are also relevant for the determination of inte- 
grating factors for higher order ODEs. As well known, any first order ODE 

N(x,u) h M(x,u) = 



11 



admits an integrating factor, i.e. a scalar function p(x, u) such that 

p [N(x,u)du + M(x,u)dx] = dF 

for some function F = F(x,u). 

The integrating factor of an n-th order equation 



d' 



N(x,u,u',...,u( n -V)— + M(x,u,u',...,u^- ,: 

is a scalar function p such that 

V u (p(M + Nu n )) = , 
where T> u is the variational derivative, 

V„ = — = D r . — + Dl 



8u du du x x du xx 

The problem of integrating factors for higher order equations 5 was tackled 
by Muriel and Romero using symmetries and A-symmctrics [49] . 

In particular, they were able to identify, given an integrating factor, a A- 
symmetry associated to it; and conversely, given a A-symmetry, a corresponding 
integrating factor. As any ODE admits A-symmetries, this implies in particular 
that - similarly to what happens for first order ODEs, as stated by the very 
classical result by Clairaut - any differential equation of order n, 

p± = /(zy™- 1 )) 

dx n JK ' 

with / analytic in some open subset of J n_1 M, admits an integrating factor 
p(x,u^ ) with some k < n. 6 



2.7 Lambda symmetries and first integrals 

As mentioned above, A-symmetries can also be considered for systems of ODEs, 
and in particular for dynamical systems. In this case, one is specially interested 
in first integrals. The relation between these and A-symmetries was considered 
by Zhang and Li [63]. 
A dynamical system 

< = f a (x,u) ; (25) 

this is represented by the dynamical vector field F = f a (x,u)(d/du a ). With 
C the Lie derivative, a vector field Z — f a (x,u)(d/du a ) is a A-Liouville vector 
field (or simply a Liouville vector field when A = 0) if 

(J^ + Cf + X^ Z + (divF) Z = . 

5 For the relation between ordinary symmetries and integrating factors, see e.g. [38]. 

6 After the first version of this paper was submitted, Muriel and Romero published another 
paper, in which they study in depth the interrelations between A-symmetries, integrating 
factors and first integrals of second order equations [50]; this also includes application to 
Ermakov-Pinney equation. Here we can just urge the reader to read this paper. 



12 



If X is a A-symmetry for (25), then (d x + Cf + X)X = 0. Moreover, let Z be 
Liouvillc for (25) and the scalar function p(x, u) satisfy the equation (dp/dx) + 
div(pF) = 0; then the vector field X := (l/p)Z is a A-symmetry for (25) [63]. 

Now, let 7 be a scalar analytic function, solution to 

(d x + C F Yi = A7 . 

Under the assumption that 7Z is divergence-free, Zhang and Li propose a way, 
starting from a Liouville vector field Z and 7, and having determined the n — 1 
invariants of the vector field Z (obtained by solving the associate characteristic 
equation), to build (n — 1) first integrals for (25). Note that this associates to 
a single symmetry not one, but (n — 1) first integrals; on the other hand, the 
requirement of determining a A-Liouville vector field and its characteristic lines 
is a rather strong one; in particular, solving for the invariants of Z is in general 
not any easier than integrating (25). 

3 The geometric meaning of lambda symmetries 

The results reviewed in the previous section are rather impressive, in particular 
in that they were able to include in the symmetry theory of differential equations 
some features which appeared, beforehand, to be definitely out of its reach. 

Nevertheless, these results were reached by an essentially analytic approach, 
and thus shed little light on the geometric meaning of lambda-prolongations 
and symmetries. Also, they were limited to ODEs; one would hope that under- 
standing the geometry of lambda-prolongation would allow for their extension 
to PDEs. In this respect, two papers were quite instrumental in bridging the 
gap between analysis and geometry of A-prolongations, and between ODEs and 
PDEs as far as A-symmetries were concerned. 

3.1 The work of Pucci and Saccomandi 

In 2002, Pucci and Saccomandi [57] devoted further study to lambda-symmetries, 
and stressed a very interesting geometrical property of lambda-prolongations: 
that is, lambda-prolonged vector fields in J n M can be characterized as the 
most general vector fields in J n M which have the same characteristics as some 
standardly-prolonged vector field. 

We stress that if Y is the lambda-prolongation of a vector field X in M, 
then the characteristics of Y will not be the same as those of the jstandard 
prolongation of X, but as those of the standard prolongation X^ of a 

generally different (for A nontrivial) vector field X in M. 

This property can also be understood by recalling (5), (7) and making use 
of a general property of Lie derivatives: indeed, for a any form on J n M, 

C XY (a) = AY~_lda + d(AF_la) = XC Y (a) + dAA(y_la) . (26) 
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3.2 The work of Morando 



It was noted [29, 45] that lambda-prolongations can be given a characterization 
similar to the one discussed above for standard prolongations; that is, with h™ 
and V as above, (20) is equivalent to either one of 

[D X ,Y] _jd = A(rji)) Vf?e£, 
[D x ,Y]=\Y + h?D m + V . 

This, as remarked by Morando, also allows to provide a characterization of 
lambda-prolonged vector holds in terms of their action on the contact forms, 
analogously to (16). 

In this context, it is natural to focus on the one- form n := Xdx; note this 
is horizontal for J n M seen as a bundle over B, and obviously satisfies D/i = 0, 
with D the total exterior derivative operator. Then, Y is a lambda-prolonged 
vector field if and only if, for all ■& g £, 

£ y (i?) + (yj%6£. (27) 



4 Twisted symmetries for PDEs: mu-prolonga- 
tions and mu-symmetries 

The result given above about the geometrical characterization (27) of A-symmetries, 
immediately opens the way to extend A-symmetries to PDEs [29]. As here the 
main object will be the one-form /x, we prefer to speak of "ti-prolongations" and 
"/U-symmetries" . 

4.1 Mu-prolongations 

Consider a semi-basic one-form 

H := Xidx 1 (28) 

on (J n M,n n ,B), satisfying D\i = 0. Then we say that the vector field Y in 
J n M fi-preserves the contact structure if and only if, for all # g £, 

C Y {$) + (7J% g £ . (29) 

Note that D/i = means DiXj = DjXi for all hence locally fi = £>$ for 
some smooth real function $. 

With standard computations [29], one obtains that (29) implies the scalar 
it-prolongation formula 

*j,< = (A + A,) * 7 - uj, m (A + Ai) r . (30) 

7 Note that the second of these relations was already remarked - and used - by Muriel and 
Romero; see eq.(22) above. 
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Let Y as in (2) be the /x-prolongation of the Lie-point vector field X (1), and 
write the standard prolongation of the latter as X^ = £'<9, + 3>j<9^; note that 
*o = $0 = f- We can obviously always write VP j = $,/ + Fj, and F = 0. Then 
it can be proved [29] that the difference terms Fj satisfy the recursion relation 

Fj. t = (A + K)Fj + XiDjQ (31) 

where Q := ip — u^ 1 is the characteristic [25, 55, 59] of the vector field X. 

This shows at once that the fi-prolongation of X coincides with its standard 
prolongation on the X -invariant space I x ; indeed, I x C J n M is the subspace 
identified by DjQ = for all J of length < | J| < n. It follows that the 
standard PDE symmetry reduction method [25, 55, 59] works equally well when 
X is a /x-symmetry of A as in the case where X is a standard symmetry of A; 
sec e.g. [22, 29] for examples. 8 

4.2 Mu-symmetries for systems of PDEs 

The developments described in the previous subsection do not include the case 
of (systems of) PDEs for several dependent variables, i.e. the case with q > 1 
in our present notation. This was dealt with in two works [22, 29], to which we 
refer for details. 

In this case the relevant contact forms are 

d a j := du} - u^dx* , (32) 

and it is convenient to see them as the components of a vector- valued contact 
form dj [60]. We will denote by & the module over q-dimensional smooth matrix 
functions generated by the dj, i.e. the set of vector- valued forms which can be 
written as ry = (Rj)l$ b j with Rj : J n M Mat(g) smooth matrix functions. 

Correspondingly, the fundamental form /x will be a horizontal one-form with 
values in the Lie algebra g£(q) (the algebra of the group GL(q), consisting of 
non-singular (7-dimensional real matrices) [60] . We will thus write 

l-i = A t dx l (33) 

where Aj are smooth matrix functions satisfying additional compatibility con- 
ditions stated and discussed below. 

We will say that the vector field Y in J n M /i-preserves the vector contact 
structure 6 if, for all G 6, 

£y(i9) + (yj(A,)^ 6 ) dx i e 6; (34) 

this should be compared to standard preservation of the contract structure in 
the form (16). 

8 The concept of /^-symmetries is also generalized to an analogue of standard conditional 
and partial symmetries [15, 20], i.e. partial (conditional) /^-symmetries [22]. 
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In terms of the coefficients of Y, see (19), this is equivalent to the requirement 
that the VP" obey the vector /i-prolongation formula 

= ( V *)?*j - < m [(V^r] , (35) 

where we have introduced the (matrix) differential operators 

V 4 := IDi + A, , 

with / the q x q identity matrix. 

We note for later reference that for vertical vector fields X = Q a (d/du a ), 
(35) yields for the coefficients of the first prolongation Y = X + ip? (d / duf) , 
simply 

C = (Vi)b<2 6 = DiQ a + (Ri) a b Q b . (36) 

If Y is the /x-prolongation of the vector field X, and Y : 5a — >• TSa, we say 
that X is a ^-symmetry for A. 

If again we consider a vector field Y as in (19) which is the ^.-prolongation 
of a Lie-point vector field X, and write the standard prolongation of the latter 
as X^ = + <S> a } d J a (with *g = % = ip a ), we can write = $ a 7 + FJ, with 
Fq = 0. Then the difference terms Fj satisfy the recursion relation 

- % [A(r J )»] {D 3 Q C ) + (Ai)j [(T J ) b c (DjQ c ) + DjQ b ] , (37) 

where Q a := (p a — uf£ l , and T J are certain matrices (see ref. [29] for the explicit 
expression). This, as for the scalar case, shows that the /x-prolongation of X 
coincides with its standard prolongation on the X-invariant space Ix', hence, 
again, the standard PDE symmetry reduction method works equally well for 
/z-symmetries as for standard ones. See ref. [29] for examples. 

4.3 Compatibility condition, and gauge equivalence 

As mentioned above the form /x, see (33), is not arbitrary: it must satisfy a 
compatibility condition (this guarantees the VP} defined by (35) are uniquely 
determined), expressed by 

[Vi.Vfc] = AA fe - DkA, + [A t ,A k ] = . (38) 

It is quite interesting to remark [22] that this is nothing but the coordinate 
expression for the horizontal Maurer-Cartan equation 9 

+ \ [M,M] = . (39) 

Based on this condition, and on classical results of differential geometry [14, 
24, 58] and a theorem by Marvan [39], it follows that locally in any contractible 

9 This expresses the requirement that the standard Maurer-Cartan equation is satisfied 
modulo contact forms, i.e. d/i + (l/2)[/i, n] £ £. 
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neighborhood A C J n M, there exists ja ■ A — >• GL(q) such that (locally in A) 
/x is the Darboux derivative of 7/1- 10 

In other words, any ti-prolonged vector field is locally gauge-equivalent to a 
standard prolonged vector field [22], the gauge group being GL(q). 

It should be mentioned that when J n M is topologically nontrivial, or /j, 
presents singular points, one can have nontrivial /i-symmetries; this is shown by 
means of very concrete examples in [22] . 

Finally, we note that when we consider symmetries of a given equation A, the 
compatibility condition (39) needs to be satisfied only on Sa Q J n M. Indeed 
when n is not satisfying everywhere (39), /i-symmetries can happen to be gauge- 
equivalent to standard nonlocal symmetries of exponential form, as noted in [41] 
and also remarked in [22]; in this respect, see also sect. 4. 5 below. 



4.4 The work of Cicogna: rho-symmetries 

The approach of /it-symmetries is inherently multi-dimensional 11 ; one could thus 
hope that it would provide better results for the analysis of systems of ODEs. 

This task was undertaken by Cicogna, who focused on systems of first order 
ODEs [17]. He considered in particular the possibility of reducing such a system; 
reduction is achieved by passing to suitable (symmetry-adapted) coordinates. 

Consider a system of n equations for u a (x), a — 1, n, given by 

F a (x,u,u x ) = (a = l,...,n) (40) 

and which admits a zx-symmetry X, with /i = Ada;; let (y, w , w n ~ , z) be 
the symmetry adapted coordinates, with y the new independent variable and z 
the invariant dependent variable. Then the system (40), when written in the 
(y, w, z) variables, depend only on the (y, w) variables and on the n first order 
differential invariants ( a , obtained solving the characteristic equation 

dz _ dw v _ dz v 

Here the are matrix functions given by 

M<°> = |J(AQ) b (o = l,...,n-l), = j^(AQ)". 

We have denoted by Q the characteristic of the vector field X: if this is written 
in the form (1), we have Q a = ip a - 

10 In a different context, it should also be stressed that the compatibility condition (38) - 
or equivalently the appearance of the horizontal Maurer-Cartan equation - also hints at a 
relation between /^-symmetries and zero curvature representations for PDEs [39], and hence 
the theory of integrable systems. This aspect of /^.-symmetries has not been studied, and could 
provide interesting results. 

^Actually the one-dimensional case is a very degenerate one, first of all because the com- 
patibility condition is trivially satisfied. 
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This results gets more clear if instead of general first order systems (40) we 
consider dynamical systems, 

< = f a (x,u) . (41) 

In this case one would like that the transformation to symmetry adapted coor- 
dinates preserves the functional form of the system, i.e. to get 

w y = 9 a (y, w, z) (a = 1, n - 1) , z y = h(y, w, z) . 

In this case, assuming again (41) admits X as a ^-symmetry with fj, = Adx, we 
actually have 

(dg a /dz) = , (dh/dz) = AfW . 

If M< a ) =0 for all a, so that no f a depends on z, the system splits into an 
(n — l)-dimensional system for the w a plus a scalar ODE z y = h(y, w, z) which 
can be seen as a "reconstruction equation" whose solution allows to pass from 
solutions to the reduced system to solutions to the full one and hence to the 
original system (41). 

This result is rather similar to what is obtained for standard symmetries, 
but it should be noted that while in that case the "last equation" reduces to a 
quadrature 12 , in this case we have in general to deal with a possibly nontrivial 
equation. When A = XI, one is reduced to the situation studied by Muriel and 
Romero [45], see sect. 2. 3. 

As stressed by Cicogna, while in general A-symmetries lead to a lowering 
of the order of the equation under study (or of one of the equations in the 
system), /z-symmetries for dynamical system do not reduce the order of the 
system, but rather split it into a reduced one and a reconstruction equation. 
Thus, he suggests that this special kind of twisted symmetry is given the name 
of p-symmetries, where p stands for "reducing" (in this context, A would also 
mean "lowering") [17]. 

4.5 Twisted symmetries and nonlocal Lie-point symme- 
tries for PDEs 

The connection between twisted symmetries and standard nonlocal ones, already 
discussed in the ODEs case (see above), is also present for PDEs. 
In this case one considers nonlocal vector fields of the form 

Pi (x,u (n ^ ) dx l ^ 

with X a local vector field in M. This is a nonlocal exponential symmetry of 
A if Z is a symmetry of A and P satisfies, at least on the solution manifold Sa, 
the compatibility condition D%Pj = DjPi. 

12 If A = 0, so that we have a standard symmetry, we obviously have M' a ' = Af( n ) = 0, 
and the reconstruction equation is a quadrature. Conversely, one can prove that if M< a ) = 
M( n ) = 0, then X is a standard symmetry [17]. 
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Then, consider the form fx = Pidx 1 . If Z is a /z-symmetry of A, and DiPj = 
DjPi on 5a, then X is a (standard) nonlocal exponential symmetry for A; and 
conversely, if X as above is a (standard) nonlocal exponential symmetry for A, 
then Z is a ^-symmetry for A [22]. 

4.6 Twisted symmetries and gauging the exterior deriva- 
tive 

It was remarked by Morando [40] that ^-prolongations and symmetries can 
also be described in terms of a deformation of the derivation operations (Lie 
derivative and exterior derivative). 

In fact, one can consider a deformed differential based on a smooth function 
/ and corresponding to a "gauging" by /: 

d d ^ := e - f d(e f f3) = d/3 + d/A/3. (42) 

It follows by straightforward computation that this is a first order differential 
operator, preserving wedge product, and satisfying d d ^ o d d ^ = (so that one 
can build complexes and a cohomology on its basis) . 

Similarly, one can consider a deformed Lie derivative based on a smooth 
function /, defined on forms as 

£%0 := e-1 £ (efx) (3 = C x f3 + df A (X . (43) 

This coincides with the ordinary Lie derivative on functions, and in general 
satisfies 

CjJ (ft+ft) = C d x f Pi + C d Jf3 2 ; 

CjJ (ft Aft) - (Cjfh) A ft + p l A (/#ft) . 

This deformed Lie derivative will be defined on vector fields as 

C d x f (Y) := e-fC (eSx) Y = C X Y - (7jd/)I. (44) 

It is immediate to check [40] this entails, for all vector fields X, Y and forms j3 
on M, 

£ d J(Y^f3) = C d J(Y)^(3 + Fj(£ d /^ . 

Moreover, Cartan formulas - with the deformed exterior derivative taking 
the place of the ordinary one - hold for the deformed Lie derivative: 

£%(J3) = X^dfi + d d f(Y^/3) , 

C d J (d(3) = d d / (C d j 09)) . (45) 

As stressed by notation, the deformed derivatives actually depend on df 
rather than on /; it is thus quite natural to define them depending on a general 
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closed (but not necessarily exact) one- form [i G A 1 (M). This yields 

d^/3 := d/3 + fi A /3 ; (46) 
C x (3 := C x [3 + fiA(X_i/3) , (47) 
:= C X Y - (Y_i/j,)X. (48) 

One can check, again by explicit computation [40] that dfi — 0, d M od M = 0, and 
££d/3 = d^(£^/3) are all equivalent. 

Needless to say, these definitions can be extended from M to J k M, i.e. to 
Jet bundles of any order, in a rather obvious way. In this case denote, as earlier 
on, by £ the contact ideal (the ideal generated by the contact forms). Then, 
d/j, e £, d^(d^) G £ and C x (dp) - d^(C x (3) G £ are all equivalent. 

Consider now a vector field Y on J n M, which leaves M invariant and which 
reduces to X when restricted to M. It can be proven [40] that Y is the /j,- 
prolongation of X if and only if the deformed Lie derivative preserves £, i.e. 

c Y iff) g £ ye G £ . 

Moreover, with T> the distribution generated by the operators D i: Y is the 
/^-prolongation of X if and only if 

Ly (Z) G T> \/Z G T> . 

This formalism is also helpful when one considers variational twisted sym- 
metries; the latter will be discussed in a later section. 



5 Twisted symmetries and gauged vector fields 

The discussion so far, in particular for ODEs, mentioned at several points the 
concept of gauge transformations. This is one of the central ideas in 20 th century 
Physics [14, 52, 53], and in a way it is surprising that it has been absent so far 
in the symmetry theory of differential equations. 13 

As stated above, see sect. 4, /^-prolonged vector fields are (locally) gauge 
equivalent to standard-prolonged ones. More precisely, if Y is the yu-prolongation 
of a vector field X, then there are vector fields W and Z, gauge-equivalent via 
the same gauge transformation (acting respectively in T( J k )M and in T(M)) 
to Y and X, and such that W is the standard prolongation of Z. This is 
schematically summarized in the following diagram: 

X Z 

fi — prol prol 
I i ' 

Y w 

For these considerations, it is convenient to deal with evolutionary representa- 
tives of vector fields [55], which we will implicitly do. 



13 It is maybe worth warning the experienced reader that the gauge transformations to be 
considered here are, in general, of a more general type than standard Yang-Mills ones. 
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We recall that the gauge group L (modelled over a Lie group G) acts in the 
same way on the vector {ip 1 , ip q } of the components of the vector field X in 
M, and on the vectors {V>j, V'j} 01 components (relative to a given multi- 
index J, i.e. to partial derivatives with respect to the same array of independent 
variables) of the vector field Y in J k M. One also says that V acts via a Jet 
representation. 

It is thus quite natural to consider a framework in which gauge transforma- 
tions are also taken into account; that is - similarly to what happens in gauge 
theories of Theoretical Physics - one introduces, beside dependent and inde- 
pendent variables, also gauge variables] these keep track of the reference frame 
changes [2 7]. 14 

At difference with the approach of gauge theories, however, the tradition in 
Applied Mathematics (in general, and in the symmetry theory of differential 
equations in particular) is to deal with standard partial derivatives rather than 
with covariant derivatives. In the specific prolongations framework, this means 
that the total derivatives appearing in the prolongation formula (5) 

1>% = D^ a j 

(for evolutionary vector field the prolongation formula reduces to this) change 
when we change frame. 

Let us choose a representation for G, and a basis {A, L r } of (left-invariant) 
vector fields for the Lie algebra Q of G; then the gauge transformations can be 
written as 

S(x) = exp[a i (a;)L i ] . (49) 

Denote by {ei, ...,e p } the reference frame for the tangent space TU to the 
manifold of dependent variables, and by {fi,...,f p } another frame, related to 
the previous one by 

f a (x) = S a b (x) e b (x) . 
A field $(x) will be given by 

*{x) = u a {x) f a (x) 

with u a the components in the f frame; in the e frame, we would have 

$(a0 = u a (x) S a b (x) e b (x) := v a (x) e a (x) . 

If now we look at x-derivatives, we get (we omit to indicate the x dependence 
for the sake of notation) 

A * = < fa 

in the f frame; while in the e frame we get, using the representation (49), 

A* = [a m (Ll) a b v\ + «r (L T m )\ v b ] . (50) 

14 It has to be stressed that here one considers changes of reference frame, and not changes 
of variables. 
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That is, field derivatives do not change in the same way as the field; as well 
known this phenomenon can be eliminated by considering covariant derivatives 
instead of standard ones [14, 52, 53]. 

We will thus consider gauge variables beside standard ones; making use of 
the basis {Li, ...,L r } we can take the a m as gauge variables. These index the 
gauge transformation, see (49) above. The function a : B — > Q is identified 
with a section of a bundle (Ag,ng, B), which is an associated bundle to Pq (the 
principal bundle of fiber G over B defining the gauge action) . 

In this way, the phase bundle is augmented from (M, ir, B) to M = (M © 
Ag,ir (Bng,B). Correspondingly one would consider vector fields 

X = C(x,u,a)^- + i P a (x,u,a) 1 ^ + A m (x,u,a)-J^ ; (51) 
ox 1 ou a oa m 

we will actually consider their evolutionary representatives, 

d d 
X = Q a (x,u,a;u x ,a x ) — + P m (x,u,a;u x ,a x )-^-^ (52) 

with components 

Q a = ip a - <f , P m = A m - a™C . 

Note these imply that the Q a will not depend on a x , nor the P m on u x . 

If X can be expressed as the gauge transformed of a vector field in M, we 
say it is a reframed vector field. It is immediate to see that such vector fields 
can be characterized as those for which Q can be written in the form 

Q a (x,u,a;u x ,a x ) = [K a b (a)\ Qq(x,u,u x ) (53) 

with K(a) G G. One can choose K(0) = I; in general we will have K{a) = 
exp[a m L m ]. We then say that X is gauge equivalent to the vector field Z = 
Q%(x,u,u x )(d/du a ). _ 

When considering prolongations in JM, the prolongation operation will also 
act on the components in the a m direction. Applying the (standard) prolonga- 
tion formula on X we will obtain a vector field 

Y = ^ 4 + ^ - &Q a ) + 0,p-) ^ (54) 

(the sum is also over all multi-indices J of suitable module | J|). 

Here the notation Dj recalls that total derivatives should be computed by 
taking into account the gauge variables a m (x) as well: 

B ' = ^ + »«4 + a "'4 :=A + z " (5f " 

If now we apply (54) on vector fields of the form (53), we get - for coefficients 
of first derivative variables - in explicit form [27] 

C = (DiQ a ) + (R t ) a b Q h (56) 
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where the matrices Ri are defined as 

Ri = [Zi{K)\ K- 1 . (57) 

These matrices satisfy 

ZiRj — ZjRi + [Ri,Rj] = ; (58) 

this is nothing else than the horizontal Maurer-Cartan equation 15 for the Ri. 
This also guarantees no compatibility problem will arise when considering higher 
order prolongations. 

It is obvious that (56) will reproduce (36) if and only if Ri = Aj. This relation 
is, however, by itself meaningless as Ri = Ri(a,a x ) and Aj = Ai(x,u,u x ). On 
the other hand, it makes sense if we consider it on a section of the gauge bundle 
{J k M lPk , J k M) with fiber p- 1 {p) = T k G p . 

A section cr 7 of the gauge bundle is defined in local coordinates by 

7 m : = a m - f m {x,u) = ; 
its lift to first derivatives is hence given by 

/Qtm df m \ 
_m . „m I J i „,a J \ ~,m ( n f'm\ r, 

7, - a, - {-Q-r + u i^r) - a i ~ W ) - °- 

Therefore, the restriction of Ri = Aj to the section a makes perfect sense, 
and reads 

A, = [aT]aL m = (Dif m )L m . (59) 

It should be noted that the restriction to a gauge section a makes sense 
provided this section is itself invariant under the action of the vector field we 
are considering; this is the case provided 

[P m - (daDQX - o . 

This can always be achieved by simply requiring P' m — (d a f m )Q a ; in facts, our 
discussion constrained Q to the form (53), but did not set any constraint on P. 
Finally, let us note that a section a of the gauge bundle (J k M,pk, J k M) is by 
definition isomorphic to J k M; thus the vector field (the restriction of X 

to cr) in a uniquely defines a vector field in J k M . 

In this way we have established a well-defined relation between reframed (via 
a gauge transformation K) vector fields on J k M, restricted to a gauge section 
a, and /x-prolonged vector fields on J k M. The form p, is given by p, — Aids*, 
and the A, satisfy 

A, = [K- 1 Zi(K)] aW . (60) 
The reader is referred to [27] for further details and other results on this line. 

15 Note that as the K only depends on a, actually ZiRj = DiRj\ hence this equation is 
cquivalcntly written as DiRj — DjRi + [Ri, Rj] = 0. 



23 



6 A gauge-theoretic approach to twisted sym- 
metries 



The discussion of the previous section (and of [27]) was based on considering 
the gauge variables as new dependent variables; as seen above, this led to some 
inconsistency, which could be cured only by restricting on a given section, i.e. 
de facto forcing the gauge variables a m to depend on (x, u) in a given manner. 

This situation is of course not satisfactory, and calls for a fully coherent - 
and fully gauge-theoretic - treatment. This was proposed in a recent paper 
[28] , and we report here the main lines of the construction proposed there. We 
also refer to [28] (see the appendix there) for a discussion of the relation of this 
approach to that of [27], and also to the approach to ^-symmetries based on the 
formalism of coverings [12]. 

In the approach discussed in the previous section, gauge variables were con- 
sidered as new, auxiliary, dependent variables and treated as such in complete 
parallel to standard dependent variables (fields) up to the point where restriction 
to a section of the gauge bundle was needed. This approach is quite unnatural 
to anybody familiar with gauge theories, as gauge variables are different than 
standard fields (matter fields in usual gauge theories [14, 52, 53]) and should be 
treated accordingly. Moreover, the gauge variables control a change of reference 
frame, which is the same for derivatives of any order - in other words, we do 
not need to consider prolongations of the gauge variables beyond order one 16 . 

The key construction is still that of gauge bundles sketched in section 5; 
by this we mean both the basic gauge bundle (p, M, M) with fiber the gauge 
group T (i.e. the set of maps 7 : M — > G), and the higher order gauge bundles 
(p fe , J k M, J k M) which also have fiber T; this is an important feature. 

When one considers that M is also a bundle (n, M, B) and that jet bundles 
J k M have several fiber structures, and in particular can be seen as bundles 
both over B in which case we write (wt, J k M, B) - and over M - in which 
case we write (o^, J k M,M) - and moreover that the same considerations hold 
for J k M, the relations among all these structures are embodied in the following 
"star" diagram [28]: 

J k M -A J k M 



W B 7T 

M -A M 



(61) 



As stressed above, we have (topologically) J k M = J k M x T; that is, the jets 
only concern standard variables, and not gauge ones. Hence, the prolongation 
operation does not involve the latter. 



16 The connection form is related to derivatives of the gauge fields through the well known 
formula A = g~ 1 dg. 
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Thus, the prolongation operation leading from M to J k M should be based on 
the usual total derivative operators Di, and hence does not involve derivation 
with respect to the gauge variables. On the other hand, a vector field in M 
will have components both in the M and in the G directions, the prolongation 
operation should be applied only to the M components. 

_As usual, it will be convenient to consider the vector bundle associated to 
J k M; in concrete terms, this means passing to consider coordinates a m in the 
Lie algebra Q of T, similarly to what we made in the previous section 17 . For 
ease of notation we will keep the same notation J k M for this. 18 

A vector field in M will then be X — £ l (x, u, a)(d/dx l )+ip a (x, u, a)(d/du a )+ 
B m (x,u,a)(d/da m ). As usual in considerations involving vector fields on jet 
bundles, it will be convenient to work with evolutionary representatives [55, 59]; 
we will consistently use these. The evolutionary representative of X is 

x = ^ = Qa £r* +pm i^> ^ 

where Q a := tp a - ufC, P m := B m . 

The coordinate expression of the prolongation X^ G X{J k M) of X is 
given by the (standard) prolongation formula; we stress once again that no 
prolongation of a m components appears. For the evolutionary representative 
Y := (X^) v = X^ k) we get, with Q} = DjQ a , 

Y = Qa >lk + Pm ^- (63) 

We are specially interested in a particular class of vector fields, i.e. those for 
which 

Q a (x,u,g;u x ) = [^(g)] a b Q b (x, u; u x ) . (64) 

These will be called, for obvious reasons, gauged vector fields. 
Using local coordinates (a;, it, a), eq. (64) becomes 

Q a (x,u,a;u x ) - [K{a)] a b 9 b (x, u; u x ) . (65) 

Here K(a) is the representation of the group element g(a) = exp(a); i.e. if G 
acts in U via the representation VP, we have K(a) = v &[cxp(a)]. Note that (64) 
and (65) do not constrain in any way the components P m of the vector fields 
along the a m variables. 

Keeping in mind that the total derivative operators do not act on gauge 
variables, we obtain immediately that 

Q} = D.jQ a = [K(a)} a b D. ; 6 b . (66) 



17 We can also think of this operation as a restriction to a neighborhood of a reference section 
in the gauge bundle, in which we can use Lie algebra coordinates. 

18 A more careful discussion is provided in [28], to which wc refer for all missing details. 
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This implies that - writing Oj — DjQ - the prolongation Y of the vector field 
X, see (63), is given by 

y=[W5A +P ™_«_. (67) 

(k) 

Now, if X a = p*X and Y = pi Y are the projection of the vector fields 
X and Y to, respectively, M and J fe M, we can state formally that Y~ is the 
p-prolongation of X n for a suitable p. In fact, we have 

X a = p*X = [K{a)]\ 6 b A ; Fo = = [tf (a)]« 6 9 b 7 ^ . (68) 

Thus Xq and Y"o are the gauge transformed - via the same gauge transformation 
- of vector fields = X and = Y such that = Y is the ordinary prolongation 
of = X . 

Note this statement is only formal, as the a variables have no meaning when 
we work in M and J k M; in order to make this into a real theorem, we will need 
to "fix the gauge", as discussed below (and similarly to what has been done 
in the previous section). Note also the relation between Xo and Yq depends 
substantially on the assumption that X is a gauged vector field. 

Fixing the gauge means selecting a section 7 of the gauge bundle. We will 
denote by u/ 7 * 1 the operator of restriction from M to M 7 , and by p^ the 

restriction of the projection p : M — > M to M 7 . We also denote by lu^ : 
J k M -> M 7 fc) and by p^ ] : M 7 fe) ->• J fe Af the lift of the maps and p^' to 
maps between corresponding jet spaces of order fc. 19 

We will summarize relations and maps between relevant fiber bundles in the 
following diagram: (y) 

M ^—^ M 7 M 



J-i) 



[j k [M] (69) 

J k M ^ M 7 fe) ^> J fc M 

As anticipated, passing to consider M 7 rather than the full M, and M 7 fc ^ 
rather than the full M^ k \ corresponds in physical terms to a gauge fixing. 

Needless to say, restriction to M 7 makes sense only if this manifold is invari- 
ant under the vector field we are considering; it is easy to check that if X be in 
the form (62), then the submanifold M 1 C M identified by a m = A m (x,u) is 
invariant under X 1 if and only if 

prn = {dA m/ dua) Q a (?Q) 

As remarked above, P m was not constrained by our previous considerations, so it 
can be adjusted to the M 7 , i.e. to the gauge fixing 7, we wish to consider. More 



19 Note that while p is of course not invertible, it follows from M 7 ~ M that is invertible, 
with (p' 7 ') _1 = 7. similarly, p( 7 ) is invertible. 
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precisely, one can show [28] that: given arbitrary smooth functions Q a (x, u, a), 
and an arbitrary section 7 of the gauge bundle, there is always a vector field X v 
of the form X v = Q a d a + P m d m such that X v leaves M 7 invariant. 

It should also be noted that X 7 projects in turn to a vertical vector field 
W on M, W — Q a (x,u)(d/du a ); and conversely any such W lifts to a vertical 
vector field W 1 = X y on M 7 , W~< = Q a (x, u)[{d/du a ) + {{dA m / du a ){d / da m )]. 

We can then consider gauging and prolongation of vector fields. In particular, 
in studying twisted symmetries one is led to consider if it is possible to find 

suitable operators and such that, for a given I on M, wc have a 
commutative diagram 

,,(7) „(7) 

X ^ x 1 ^ w 



Prfc[M] 

xik ) <y>i_i x 



(71) 
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We will just give the main results, referring again the reader to [28] for 
a complete discussion. We will introduce, in order to state our result in a 
compact form, operators := owW, r^ 7 ' : M — > M; and correspondingly 
:=^o W W,tW:M( 1 )4MW. 

( k) 

Then we have that: The twisted prolongation operator P 7 is uniquely de- 
fined by the requirement that (t^ 7 ' 1 )* o (Prfc[M]) = o Moreover, P 7 fe ' 
corresponds to the /i-prolongation operator of order k with [i = [d\I/(7)]\E'(7 -1 ). 
With the local coordinates (x,u,a), this corresponds to fi = Aidx 1 where Aj = 
-R^ = -(D.K^K- 1 = K 1 (D i K~ 1 ). 

A converse of the result above is as follows: Let Y be the k-th ^.-prolongation 
of the evolutionary vector field W on M, with p 6 A 1 ( J 1 M, ip(G)) given in 
coordinates by fi — Aj(x, u, u x )dx l . Then: (i) there is a section 7 of the gauge 
bundle such that Y = V 1 (W). (ii) there is a vertical vector field X in M such 
that (71) applies, (iii) The matrix function K 7 (x,u,u x ) satisfies DiK 7 — 
-K l K 1 . 

To conclude this (long) section, let us briefly comment on the geometrical 
meaning of the main results just given. We have shown that the ^-prolongation 
operator appears if we are insisting in restricting our analysis to the phase bundle 
M (or to the sub-bundle M 7 C M seen as an image of M under the gauge map 
7 embedding it into M) rather than to the full gauge bundle M. The fact we are 
considering projections of vector fields in M 7 C M and C to vector 

fields in M and makes that the relation between basic vector fields and 

prolonged ones is not the natural one, described by the prolongation operator, 
but is the "twisted" one described by the /x-prolongation operator. See also the 
discussion in the Appendix to [28]. 
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7 Twisted symmetries and variational problems 



As well known, symmetry analysis - and symmetry reduction - are specially 
powerful in dealing with variational problems; the key result relating standard 
symmetry and reduction is in this framework the classical Noethcr theorem [54] , 
see also [35, 55]. 

It is thus not surprising that the application of twisted symmetries to varia- 
tional problems is specially fruitful and - in this author's opinion - also specially 
fascinating. 

Unfortunately, this paper is already way too long - which also accounts for 
the lack of detailed examples in it - so that we will just very sketchily mention 
the sources of the main results obtained in this direction, leaving to the reader 
to look directly the original papers 20 . 

The first study in this direction was conducted by Muriel, Romero and 
Olvcr [51]; they considered variational problems defined by a Lagrangian (of 
arbitrary order) in one dependent and one independent variable, and studied 
both how A-prolongations allow to construct new methods for the reduction 
of Euler- Lagrange equations and a version of Noether's theorem adapted to A- 
symmetries. In this case one focuses on "variational A symmetries" - which are 
the twisted symmetries analogue of standard variational symmetries [55] - and 
obtains partial conservation laws. With L the Lagrangian, X a vector field and 
Y its A-prolongation, X is a variational A-symmetry of L if 

Y(L) + L(D + \)i = (D + \)B (72) 

for some function B(x, u, u x , ...). 

If the n-th order Lagrangian L admits a variational A-symmetry, then there 
exists a Lagrangian L of order (n — 1) such that a (2n — l)-paramctcr family 
of solutions to the variational problem described by L can be found from the 
solution to the variational problem described by L by solving an auxiliary first 
order equation. The result extends to generalized variational A-symmetries. 

Moreover, denote by E\L\ the Euler-Lagrange equation corresponding to L. 
If L admits a variational A-symmetry X, and Q is the characteristic of X, then 
there is some P(x, u, u x , ...) such that 

Q E[L] = (D + \)P . (73) 

It turns out [51] that in this case X is a A-symmetry for the equation P[u] = 
0; the reduced equation obtained from P = using the symmetry X is, up to 
multipliers, the reduced equation E[L] for the variational problem (see above). 

Finally, if one restricts on the solutions to the Euler-Lagrange equation 
E[L] = 0, then (1/P)X is a variational symmetry of the variational problem 21 . 
This allows to associate a partial conservation law to such symmetries. We refer 
again the reader to [51] for details. 

20 We hope to be able at a later time to also review developments in this direction. 
21 When this happens, i.e. we need the restriction to solutions of the Euler-Lagrange equa- 
tion to have a variational symmetry, one speaks of "pseudo- variational symmetries" [55]. 
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It should be stressed again that the Muricl-Romcro-Olvcr approach is not 
restricted to first order Lagrangians, and they consider explicitly some higher 
order example. On the other hand, as implied by consideration of A-symmetries, 
they consider Lagrangian problems with only one independent variable. 

In her work on deformation of the Lie derivative [40], Morando considered 
variational A- and ^-symmetries. She started from a geometrical characteriza- 
tion - in terms of the action of the deformed Lie derivative (see section 5 above) 
on the Poincare-Cartan form O - of variational A-symmetries, which reads 

££(6) e £ ; (74) 

here Y is the A-prolongation of the vector field X on M, and £ is the contact 
ideal (again, see above). 

This characterization is immediately generalized, and written exactly in the 
same form, to yu-symmetries and hence field theory. 

To a divergence variational A-symmetry AT of a first order regular Lagrangian 
L is associated the "A-conservation law" 

D x (X_l6 - R) + A(X_|6 - R) , (75) 

where R is a suitable smooth function. 

In the case of field theory, it is convenient to introduce a form 

p = R l fl t + a , 

where a G £ and Cti = di _l Ct, with O = dx 1 A ... A dx p the volume form on the 
base manifold B. In this case one has a "^(-conservation law"; this reads 

Di _l d^pf _l6 - p) = (76) 

in compact geometrical notation; in coordinates, and using freely the standard 
notation introduced earlier on, it yields 

(A + Ai) (H - uie) + eL - Ri ) = ° • (77) 

In a different work [21], Cicogna and Gaeta also looked at /i-symmetries for 
variational problems; they also obtained that to a /i-symmctry of a Lagrangian 
is associated a deformed conservation law as above. 

On the other hand, they also looked at twisted symmetries of variational 
problems from the point of view of gauged vector fields. In this framework, it is 
natural to wonder how the Euler-Lagrange equation themselves are transformed 
by a change of reference frame. One thus obtains "/i-Euler-Lagrange equations" , 
and this are exactly invariant under (the ^(-prolongation of) a vector field which 
is a /i-symmetry for the underlying Lagrangian; correspondingly, a /i-symmetry 
for the Lagrangian yields an exact conservation law for the associated /x-Euler- 
Lagrange equations [21]. 

All these works were conducted in the Lagrangian framework; in a very 
recent study, Cicogna considered the Hamiltonian counterpart of these results. 
We will not discuss at all this aspect, referring the reader to his paper [18]. 
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